Coherent Control for a Two-level System Coupled to Phonons 
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The interband polarizations induced by two phase-locked 
pulses in a semiconductor show strong interference effects de- 
pending on the time n separating the pulses. The four-wave 
mixing signal diffracted from a third pulse delayed by r is 
coherently controlled by tuning t\ . The four- wave mixing re- 
sponse is evaluated exactly for a two-level system coupled to 
a single LO phonon. In the weak coupling regime it shows 
oscillations with the phonon frequency which turn into sharp 
peaks at multiples of the phonon period for a larger coupling 
strength. Destructive interferences between the two phase- 
locked pulses produce a splitting of the phonon peaks into a 
doublet. For fixed r but varying n the signal shows rapid 
oscillations at the interband-transition frequency, whose am- 
plitude exhibits bursts at multiples of the phonon period. 



Ultrafast spectroscopy on semiconductors gives a di- 
rect insight into the loss of coherence between optically 
excited carriers. The dephasing takes place through dif- 
ferent mechanisms such as emission of phonons, Coulomb 
interaction or scattering by impurities. In the first hun- 
dred femtoseconds after an optical excitation few scat- 
tering processes have occurred and the electronic wave- 
function may be described by a superposition of different 
states with well defined phase relations. The resulting 
quantum beats are a signature of phase coherence. 

Quantum beats with the longitudinal optical (LO) 
phonon frequency have been observed in four-wave mix- 
ing (FWM) experiments on GaAs which was excited res- 
onantly at the excitonic energycl. The beatings are due 
to interference between states with different number of 
virtual phonons. Indeed electrons close to the band edge 
which do not have enough energy to emit any real LO 
phonon, form a polaronic state. 

The phase coherence of the excitation may also be used 
to coherently control the optical response of the sample 
via a pair of phase-locked pulses which induce interfer- 
ing polarizations in the sampleB. Recently this coherent 
control was applied to the quantum beats in the FWM 
signal]. Depending on the time delay t\ between phase- 
locked pulses the amplitude of the quantum beats was 
modulated, thus experimentally demonstrating that the 
phase coherence is not lost in the first few scattering pro- 
cesses. 



The present work computes the FWM response in third 
order in the field for a simple two-level system treating 
exactly the coupling to a single phonon mode. The cal- 
culations closely follow Ref. which studied the linear 
response in connection to phonon broadening of impurity 
spectra in semiconductorsB. Recently the FWM response 
in the same model was obtained for (5-function pulses by 
solving the infinite hierarchy of kinetic equationsO. 

We focus here on the regime of strong coupling between 
electrons and phonons and discuss the recent coherent- 
control experiment on the polar material ZnSe, where 
peaks separated by approximately half the phonon pe- 
riod were observed and interpreted as signatures of two- 
phonon processes! In the two-level system, however, we 
attribute these features to a splitting of the phonon peaks 
into doublets due to destructive interference between the 
signals generated by the two phase-locked pulses. We 
further predict non-perturbative bursts of the FWM sig- 
nal as a function of delay between phase-locked pulses at 
multiples of the phonon period T p h- 

A two-level system | i), i = 0, 1 with energies eo < t\ 
couples to the position a + of a harmonic oscillator of 
frequency fl with a level-dependent coupling strength gi. 
In addition the classical light field E(t) induces transi- 
tions between the two levels: 

H = £i=o,i ( e * + 9i(a + al)) | +n a ^a 

-f,E(t)(\0)(l\+H.c). (I) 

This model describes the coupling of electrons to LO 
phonons of frequency fl — 2%/Tph using a single oscilla- 
tor. The two levels may be viewed as states in the va- 
lence and conduction bands which are resonantly excited 
by light. Within this interpretation the coupling con- 
stants go and g\ should be taken equal since the Frohlich 
interaction does not depend on the effective mass of the 
electrons i.e. on the bandtl Another picture associates 
the level | 0) to the full valence and empty conduction 
band while | f) corresponds to an excitonic state which 
has an effective couplingd gi to the phonons differing from 
the ground-state value go. This interpretation is partic- 
ularly meaningful for the FWM experiments where the 
excitation is resonant with the excitonic level. 

The model considers only scattering of the electron 
within the same level and hence cannot describe relax- 
ation effects. This approximation which neglects the re- 
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coil of the electron is particularly relevant to the strong 
coupling regime since it corresponds to the small-polaron 



limitu£3. 

Four-wave mixing. A first pulse E\(t) creates an in- 
tcrband polarization in the sample while a second pulse 
E2 (t — t) delayed by a time r is diffracted along the di- 
rection 2q2 — qi where qi and qi are the propagation di- 
rections for the first and second pulses, respectively (see 
inset of Fig.Q). The FWM signal F{t,r) atJJme t and 
delay r is related to the third order response^ x^ 3 -*: 
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E 2 {ti - T)E 2 {t 3 - t) X {3) (t-t 1 ,t-t 2 ,t-t 3 ) (2) 

The electron occupies the lowest level | 0} before the 
pulses arrive, and the oscillator is in thermal equilibrium 
at an inverse temperature (3 — 1/ksT. The pulses cause 
transitions between electronic levels which affect the evo- 
lution of the oscillator via the level dependent Hamilto- 
nian Hi = fla^a + + gi(a + a^). The response function 
involves a thermal average of evolution operators and has 
a main contribution at positive time delays and a second 
term for t\ > t 2 > t 3 i.e. when the two pulses overlap: 

X {3) (ti,h,t 3 ) = e(ti)9(t 2 -ii) x 

^Q(f 3 ^ e iH i( t 2-t 1 ) e iHot le -iH 1 t 3e iH (t 3 -t 2 )^ _|_ 

e(t 3 ~t2)(e~ zHltl e lHo( - t2 - f ^e lHl{t2 - t3) e lHof3 )y (3) 

The evaluation of the terms in brackets follows ex- 
actly Ref. ||. The Hamiltonian Mi is diagonalized 
by the Lang-Firsov transformation^ S — g(a' — a): 
exp(S')iJi exp(— S) = Hq + e where e = t\ — eo + 
(<7q — 9i)/Q is the renormalized transition energy and 
g = (gi — go) I SI a dimensionless coupling constant. Since 
only the difference g\ — go enters the transformation, it is 
essential that e.g. excitonic effects renormalize the cou- 
pling strengths from the non-interacting value g\ = go. 

The product of evolution operators in Eq. (^) is evalu- 
ated using the previous canonical transformation and the 
following relation between the time-dependent operators 
S{t) = exp(iH t)S exp(-iH Q t): 



3 S(ti) -S(t 2 ) _ ^(t^-Sfa) e -ig 2 sm(Q(ti-t 2 )) 



(4) 



We do not reproduce here the lengthy expression for 
the response function but instead discuss the FWM signal 
F(t,r) for ^-function pulses. The signal oscillates at the 
renormalized energy e with a time delay 2t typical of a 
photon echo. The polaronic nature of the state shows up 
in the exponential dependence on the coupling constant, 
and the temperature enters via the Bose function Nb = 
f/(exp(/30) - 1): 

F{t,r) = -ie(r)Q(t - T)e" l£{t - 2T) x 

e -g 2 ((A r B + l/2)|2e is " !T -e is "-l| 2 -2isin(ar)+isin(Sli)) (5) 



The FWM signal which depends simply on fir exhibits 
quantum beats at harmonics of the bare LO frequence, 
in contrast to the frequency shift seen experimentallytl. 
While the intraband scattering is missing in the present 
exactly solvable model, the simulations of Ref. [I] for a 
full two-band model do reproduce the correct frequency 
within approximate quantum kinetic equations. More- 
over the FWM signal in Eq. (H) has evenly spaced phonon 
replica in the frequency domain, which are not observed 
in experiments. 

Coherent control. Two phase-locked pulses of equal 
amplitude propagate in direction q\ and are separated by 
a time delayll3 t\ as shown in the inset of Fig. [I]. Since 
the FWM is linear in the field Ei(t), the total signal 
is simply F(t, r) + F{t + t\, t + ti) which shows strong 
interference effects depending on t\ through the relative 
phase of the induced polarizations. 

The time-integrated FWM /(r,n) = \ F{t,r) + 
F(t + 7*1, r + t%) 1 2 dt can be calculated analytically for 
(5-function pulses, E\(t) = 5(t) + 5(t + Ti). A phenomeno- 
logical lifetime T which accounts for the damping of the 
polarization is introduced as an imaginary part of e. For 
simplicity we give here only the result at zero tempera- 
ture (Nb = 0) and for T <§; J7, which involves the modi- 
fied Bessel function Iq: 

J(r, n ) =G( T ,0)+G(T + Ti,0) + 2G(T + y,ri). (6) 
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We first discuss the contribution of a single pulse when 
the integrated FWM equals G(r, 0) which decays expo- 
nentially with time delay and is modulated by terms 
oscillating with multiples of £1. In the weak coupling 
regime, only the first harmonics cos(rir) contributes to 
the modulation as given explicitely in Ref. ^[ Figure |l| 
shows the integrated FWM response I(t, t\ = 0) calcu- 
lated with the exact x^ 3 ' for pulses of finite width which 
are taken as sech 2 -shaped pulses resonant jvith the two 
levels and of 13fs duration. The parameters! for ZnSe are 
fl = 32meV, e = 2.8eV, T = 77K and 1/T = 300fs. The 
weak modulation of the signal for a small coupling g = 
0.2 relevant to GaAs, evolves into distinct peaks at multi- 
ples of the phonon period for a larger coupling g = 0.8 as 
in ZnSe. In the strong coupling regime indeed, the many 
harmonics sum up into sharp peaks of width proportional 
to 1/^/g: G(t,0) ~ exp(-Fr - 16g 2 sin 4 (f7r/2)). 

Two phase-locked pulses in direction q 1 produce inter- 
ference effects depending on the delay n . When the two 
levels are decoupled from the phonons (g = 0), the FWM 
signal in Eq. (^) oscillates as sin 2 (eT!/2), and is thus sup- 
pressed when the induced polarizations are exactly out 
of phase, i.e. en = (2n + 1)tt. 



2 



In the strong coupling regime, the first two terms on 
the right-hand side of Eq. correspond to the separate 
responses of each pulse which are peaked at r = nT p h 
and nTph — Tl , respectively. The third term describes the 
interference between polarizations and has its maximum 
at r = nTph — Ti/2 in between the two other peaks. De- 
pending on its phase whose main contribution is eri , the 
interferences will be destructive or constructive. 

Figure ^ shows the integrated FWM signal as t\ is 
tuned within a narrow range around n = 4l7r/e = 
0.23T ph . At time n = 0.227T ph the different contribu- 
tions in Eq. (||) add up constructively to form a single 
broad peak at r = T p h — T\/2. At t\ = 0.231T p ^, how- 
ever, the destructive interference splits the peak into two 
bumps separated by a time r ~ 0AT p h- The separation 
exceeds T\ because the interference suppresses the sig- 
nal in the whole range of r where the peaks G(r, 0) and 
G(t + n, 0) have a significant overlap. 

The splitting of the phonon peaks was interpreted in 
Ref.0 as a doubling of frequency due to two-phonon pro- 
cesses since the peaks were separated by approximately 
Tph/ 2 for the values of n investigated. Within the two- 
level model, however, the additional peak that appears 
due to destructive interference, clearly moves with vary- 
ing ti and cannot be interpreted as a frequency doubling. 
Figure || shows the integrated FWM signal for various T\ 
tuned as to produce well separated doublet peaks, P~ 
and P + , for time delays close to eri = 31,41,51 and 
6l7r. The P~ peak follows the dispersion T p h — t± as il- 
lustrated by the inset which shows the positions of the 
maxima r max as a function of t\ for coupling strengths 
g = 0.8 and 1.5. While for g = 0.8 the rather broad peaks 
are splitted by more than t%, for g = 1.5 their positions 
approach the expected values of T max = T p h — t\ and T v h 
because the peaks have a reduced overlap. 

The FWM signal for fixed r but varying t\ shows 
bursts at multiples of the phonon period which are a clear 
signature of the strong coupling regime. Figure ^ shows 
the intensity of the first phonon peak I(r — T^n) 
rapidly oscillating with the interband frequency e. The 
envelopes J±(ti) of the oscillations shown as thick lines, 
vary exponentially with x — sin 2 (f2Ti/2) for large g: 

I±{ti) oc 1 + cxp (-2IYi - 2.g 2 (l + Ax - y/1 + 8xj) ± 
2exp(-rTi - 2.g 2 (l + 2x - (1 + 2x + \[\ + 8.x)/2) 1/2 ) 

For large g the envelopes exhibit bursts at multiples of 
the phonon period. The bursts are a non-perturbative 
effect involving all the many harmonics that contribute 
to the FWM signal. For an arbitrary time delay t\ the 
harmonics which are out of phase destroy the phase co- 
herence induced by the pulses, and the oscillations are 
suppressed exponentially rapidly with t\. For resonant 
times Ti however, the different harmonics add up coher- 
ently and the phase sensitivity is recovered. 

In conclusion we have studied the FWM response of a 
two-level system coupled to a single phonon mode focus- 



ing on the strong coupling regime. Analytical calcula- 
tions of the integrated FWM signal provided new insight 
into the interplay between coherent control and dephas- 
ing through scattering processes: (i) The FWM signal 
has sharp peaks at multiples of the phonon period whose 
width decrease as l/^/g with increasing coupling, (ii) 
Two phase-locked pulses delayed by ti give rise to two 
sets of peaks which add up coherently into a single set 
of bump at nT p h — t±/2 for en a multiple of tt. When 
the two signals are out of phase, however, the peaks split 
into doublets at approximately nT p h — t± and nT p h- This 
picture contrasts with the interpretation of the coherent 
control experiments in terms of a frequency doubling, 
(iii) The FWM signal for fixed t exhibits bursts of coher- 
ence at resonant values of Ti = nT p h, which are a clear 
signature of the non-perturbative regime and might be 
observed in experiments by measuring the intensity of 
the first phonon resonance as a function of t\ . 
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FIG. 1. Normalized integrated FWM 

signal I{t,t\ = 0)/J(0,0) as a function of delay r in units 
of the phonon period T ph for a single pulse in direction qi, i.e. 
T\ — 0. The weak modulation of the signal for g = 0.2 turn 
into distinct peaks at multiples of the phonon period with in- 
creasing coupling. The inset shows the setup for the coherent 
control where two phase-locked pulses delayed by n propa- 
gate in direction qi. A third pulse in direction §2 delayed by 
t is diffracted along 2q2 — qi . 
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FIG. 2. Integrated FWM signal J(r, n) as a function of 
delay r for different n around en ~ 417T. Constructive inter- 
ferences between phase-locked pulses produce a single peak at 
t = Tph — n /2 which splits into a doublet when n is tuned 
to achieve destructive interferences. 
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FIG. 3. Integrated FWM signal I(t,ti) for g = 0.8 as a 
function of time delay r for different times n . The inset shows 
the position of the maxima T max of the two peaks P~ and 
P + as a function of n for g — 0.8 and 1.5. With increasing 
coupling, the dispersion of the P~ peak approaches T ph — r\ 
plotted as the solid line. 
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FIG. 4. Normalized intensity of the first phonon peak 
I(t — Tph,Ti) I I(t = T p h,0) as a function of delay time n 
between phase- locked pulses. The rapid oscillations of fre- 
quency e which are first suppressed with increasing n, show 
bursts at later times when n is a multiple of the phonon pe- 
riod Tph- The envelopes of the oscillations are well described 
by Eq.(8) valid for large g (thick lines). For better display of 
the interference patterns, the frequency e has been reduced 
by a factor of 4 from its value for ZnSe. 
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